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DISTRIBUTION  OF  DROPS  ACCORDING  TO  SIZE  IN  ATOMIZING  LIQUIDS 
BY  KEANS  OF  CENTRIFUGAL  ATOMIZERS 
by 

L.  G,  Golovkov 

The  article  presents  further  support  of  the  function  of  distri¬ 
bution  of  drops  according  to  size  proposed  in  the  work  (5),  explains 
the  method  of  parameter  derivation  for  this  function  and  establishes 
the  numerical  value  of  its  parameters  for  centrifugal  atomizers. 


The  most  fruitful  concept  for  the  establishing  of  analytic  dependences  between 
parameters  which  determine  the  process  of  atomization  and  sizes  of  the  spectrum, 
in  oiir  view  proves  to  be  the  concept  expressed  for  the  first  time  by  Mugele  and 
Evans  (4)  and  developed  afterwards  by  Tresh  (5).  Its  basis  is  composed  by  two  pro¬ 
positions:  a)  in  the  spectra  of  the  atomized  liquid  there  are  always  contained  drops 
of  the  greatest  size  b)  dn  ps  of  all  the  remaining  sizes  occur  as  a  result 

of  statistical  conformity  to  rule. 

If 

Tresh  explains  the  statistical  character  of  the  functions  of  distribution  jf 
dr^ps  according  to  size  by  the  fact  thf t  at  the  original  moment  of  scattering  of 
the  Jet  there  occurs  a  disordered  .  ..change  of  energies  among  the  separate  particles 
of  liquid  through  the  irerging  of  ir iividual  drops,  their  secondazy  atomisation,  re¬ 
peated  merging,  etc.  This  chaotic  process  also  results  in  the  formation  in  t'i.  .  -vr- 
rum  of  the  atomized  liquid  of  only  a  small  number  of  very  big  and  very  small  drops. 

Hej.*e  the  author  proposes  the  folic  /ing  dependences  for  computing  the  spectra 


of  the  atomized  liquid:  the  equation  for  determining  D. 


max 
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The  expression  (2)  Is  a  well-founded  aemiempirical  function  of  the  distribut¬ 
ion  of  the  drops  by  sizes  and  consequently,  it  is  favorably  distinguished  from  the 
other  purely  empirical  equations  (for  example  (1-3,  6)  and  of  other  authors).  How¬ 
ever,  the  equation  (l)  and  the  numerical  value  of  the  parameter^  (equal  to  0.35) 

n 

were  established  by  Tresh  for  developing  the  experiments  In  pneumatic  atomization 

of  liquids  and  therefore  proved  to  be  not  suitable  for  centrifugal  atomizers. 

In  this  article  there  is  analyzed  a  supplementary  basis  of  the  function  of  dis¬ 
tribution  proposed  by  Tresh  (2),  methods  are  given  for  the  determining  of  the  para¬ 
meter^  of  this  function,  there  is  established  its  numerical  value  for  centrifugal 
atomizers,  and  there  is  also  refined  the  generalized  dependence  for  determining 
The  function  (2)  was  established  by  Treshy  the  methods  of  statistical  thermo¬ 
dynamics  with  the  postulation:  l)  of  constancy  of  the  overall  volume  of  the  drops; 
2)  of  constancy  of  the  surface  of  the  drop  mixture  in  the  period  of  chaotic  exchange 
of  energies  among  the  particles  of  the  liquid.  If  the  first  assumption  causes  no 
particular  uncertainties,  then  the  second  does  not  prove  to  be  manifest.  Below 
there  is  given  the  proof  of  TrOsh'a  second  postulate. 

The  known  peculiarities  of  the  conditions  of  the  existence  of  molecules  on  the 
boundary  separating  liquid  and  gaseous  phases  are  the  reason  that  the  increase  of 
the  surface  of  the  separation  requires  the  expenditure  of  work.  The  amount  of  this 
work  referred  to  unit  of  area  of  surface  represents  the  surface  tension  or  the 
specific  isobar  potential  Z  (7). 

In  the  formation  of  the  surface  by  the  area  ^  f  the  change  in  the  isobaric  po¬ 
tential  amoxints  to 

AZ  =aA/. 

(3) 

On  the  other  hand  it  is  known  that  for  any  process  (  3  Z/^T)p  =  -S.  Consequently, 

the  change  in  the  entropy  in  this  process  is  equal  to  /d  S  ^  Z)/  3  T)^,  or 

with  the  taking  into  consideration  of  the  relationship  (3) 
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AS 


dT 


By  disregarding  the  change  in  the  surface  corjiected  with  the  change  in  the  tempera¬ 
ture  (in  view  of  the  insignificant  coefficient  of  the  volumetric  expansion  of  liquids) 
we  shall  get 


Now  let  us  assume  that  in  the  process  under  consideration  of  the  chaotic  ex¬ 
change  of  energies  oetween  the  particles  of  the  liquid  their  total  surface  Increases 
oy  (  f  ^0).  Then  according  to  the  relationship  (3)  end  (4)^  Z  0  and  ^  S 

^0  inasmuch  as  (  *)p  a  .ways  is  negative  (with  increased  temperature  the 

Surface  tension  always  becomes  less). 

Consequently,  with  f  ^  0  the  indicated  pzocess  should  be  accompanied  simul¬ 
taneously  both  by  increase  in  the  isobarlc  potential  Z  and  oy  increase  in  the  antropy 
S.  According  to  the  second  law  of  thermodynamics  the  simultaneous  Increase  of  thes»,‘ 
two  values  in  a  spontaneous  procesf  is  not  possible.  Therefore  the  value  f  can¬ 
not  be  positive. 

The  proposition  of  the  reduction  in  the  total  surface  of  the  drops  {  ^  O) 

also  leads  to  a  contradiction  with  the  second  law  of  thermodynamics  {  ^  C  uni 

.multaneously  ^  S^C).  The  contradiction  disappears  when  f  =  0,  i.  e., 
w-.cn  the  surface  of  the  droo  mixture  ii  the  period  of  its  "disordered"  state  is  un¬ 
changed,  which  was  what  was  to  be  proved. 

In  the  solution  of  each  statistical  problem  after  establishing  the  raathtform  of 
the  function  of  the  distribution,  it  is  necessary  to  develop  methods  of  determining 
the  numerical  valuea  of  the  parameters  of  this  function  on  the  basis  of  the  exlst- 
In-  results  of  observations,  i.  e.,  to  derive  its  parameters.  The  parameters  of 
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the  function  (2)  are  and  . 

In  the  process Ing  of  the  experimental  data  on  the  atomization  of  liquids  by 
centrifugal  atomizers  with  the  aid  of  the  theory  of  dimensionalities  there  was  ob¬ 
tained  the  following  generalized  equation  for  determining 

with  regard  to  the  parameter^ as  depends  on  two  basic  methods  of  primary  pro¬ 
cessing  of  experimental  data  on  th  atomizing  of  liquids,  there  are  proposed  below 
two  operations  for  its  determination. 

For  the  case  where  to  measurement  In  the  sample  there  Is  subjected  each  of  Its 
eloments  (every  drop)  the  determination  of^is  done  by  the  method  of  the  maximum 
plausibility  (8)  which  gives  the  most  nearly  “precise"  value  for  the  parameter  , 

In  determining  in  the  sample  the  number  of  dr  ps  which  fall  Into  different 
Intervals  as  to  size  for  finding  the  numerical  values  of ^one  uses  the  method  of 


moments • 


Method  of  Maximum  Similitude.  The  relative  frequency  of  the  drops  dn/nQ  in  the 
formula  (2)  can  be  treated  as  the  probability  of  the  occurrence  of  drops  of  dimen¬ 
sions  of  from  ^  to  5  ^  b  t  and  therefore  the  function  (2)  is  the  probability  den¬ 

sity  of  the  random  value  By  accomplishing  on  the  equation  (2)  the  functional 
conversion  D  =  ^  ,  we  shall  get  the  function  of  the  density  of  the  distribution 

of  directly  in  accordance  with  the  diameter 


9(D)- 


dn  1  / 

/  1  pvn  1  - 

P  Dn,,, 

\dD 

tyiD  “  Z3*  V' 

D  ) 

] 

•  0 

D 

Jd* 

(6) 


or  after  integration 


9(D) 


eXP  (—  P  Dni>»/D)  _  , 
+P)exp(-1>) 


(7) 
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Wfc  propose  now  that  we  have  a  sample  of  the  volume  0,^  consisting  of  n  indepe:  '’  ■'t 
values  for  the  dimensions  of  the  drops 

DiJ  D2J  D3 

The  probability  of  the  occurrence  of  a  given  sample 

P{0„-,  ?)*=?(£»,;  p)<p(D,:  ?)...<?(£>„;  ?).  (8) 

representing  a  function  of  the  fixed  values  of  the  dimensions  of  the  drops  D2, 

Dj,  *.•  ,  Djj  and  the  variable  parameter^ is  also  called  the  function  of  similitude. 

It  is  apparent  that  if  the  selected  value  of^ls  such  that  P(0q;  ^  )  is  very 
small,  then  this  event  "practically"  should  not  occur  in  a  single  experiment.  Con¬ 
sequently,  such  an  hypotheses  about  the  value  oC ^  should  be  rejected  as  not  pro¬ 
bable. 

The  method  maximum  probability  applicable  to  the  given  case  consists  in  this 
that  the  "best"  value  for  is  determined  in  a  way  as  the  variable  argument  of  the 
function  (8)  at  which  it  reaches  the  maximum  value. 


In  other  words,  the  "best"  value  for ^ we  get  from  the  condition 

P  (Oj-.;  ^ )  =  max 

or,  what  is  the  same,  from  the  condition 

P  ~  ^(^1»  ^2*  ***  ’  ^n’/^  ^  ~ 

inasmuch  ns  P(0jj;  ^  )  ^0. 

The  function  of  probability  for  the  distribution  (?)  in  accordance  with  the  equation 
(8)  has  the  form 


P(0n\?)' 


D\D^...Dl(\  -i  ?rexp(-np) 


and  therefore 


«  2/J  In  P  -4-  In  —  t* 


on  V_L_y 

-  P  ,  ry  /  J 


D, 
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The  value  sought  Is  determined  by  the  condition  of  the  extremum  ^  ^dB  =0, 


hence 


P  l+P 


ms.  t 

Om,,  ^  _}_ 

n  /r 


- 1. 


After  designating  the  right  side  of  the  expression  obtained  by  ]2  we  shall  have 

P=*(l— 6  +  /(!  -6)*  +  86)/2i».  (9) 

Thus,  if  one  measures  in  the  sample  frum  the  drop  mixture  the  dimensions  of 
each  individual  drop,  then  the  finding  of  the  parameter  is  reduced  to  very  simple 


operations,  the  computation  of  the  value 


n  ^  D. 

<=i  ‘ 


-1 


and  the  use  of  the  formula  (9)« 

In  this  situation  should  be  determined  in  accordance  with  the  criterional 
equation  (5)* 

Method  of  Moments.  The  essence  of  this  method  of  deriving  parameters  consists 
In  this  that  according  to  the  Khlnchln  theorem  with  a  great  volume  of  sampling  0^ 
from  the  general  aggregate  of  objects  any  theoretical  moment  of  the  order  i  of  ran¬ 
dom  magnitude  is  approximately  equal  to  the  corresponding  empirical  moment.  In 
this  situation  the  greater  the  volume  of  the  sampling  the  more  nearly  accurate  this 
apporxlmate  equality  will  be.  For  our  case  the  most  convenient  is  to  use  the  first 
initial  moments  of  the  random  value  D,  or  what  is  the  same,  the  theoretical  Mj,(D) 
and  the  emperlcal  M^  (D)— the  mathematical  expectancy  of  the  random  value  D. 

Inasmuch  as  the  random  value  of  0  (diameter  of  the  drop)  changes  from  0  to  D, 
and  its  probability  density  is  given  by  the  equation  (6),  therefore 


max  j 


dP 

CP 
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or  after  integration 


With  the  aecond  method  of  processing  the  primary  data  cn  the  atomization  the 
results  of  the  experiment  present  themselves  thus:  in  the  intervals  of  grouping 
0— the  numbers  of  drops  respectively,  are  nj^,  Ug,  n^,. 
For  the  composition  of  (D)  let  us  have  recourse  to  the  following  operation.  Let 
us  designate  by  the  average  of  the  ith  interval  of  grouping,  ^  the  width  of 
this  interval,  and  by  ^  l(*^l)  the  relative  frequency  of  drops  of  the  size  Dj  ^  0.5 
^  Dj.  Then  f  i(D^)  =  ni/nQ  ^  and 


imX 


in) 


On  the  basis  of  Idiinchin's  theorem  Mj.(D);55:?M^  (D),  hence 
by  designating 


yD,^t(D,)^D, 


we  finally  get 


'?  =  a/(I  —a). 


(1- 


(.3) 


Consequently,  the  determination  of^  in  this  case  is  reduced  to  the  computafL'  a  .f 

the  value  and  the  application  of  the  formula  (13).  In  this  situation  D 

max 

should  be  computed  in  accordance  with  the  equation  (5). 

At  our  disposal  there  were  experimental  data  on  the  spectra  of  the  atomi  aation 
of  liquids  in  the  form  of  the  above  indicated  tables.  Therefore  the  numerical  val¬ 
ue  for  the  pararaetery^ for  the  centrifugal  atomizers  was  determined  by  the  method 
of  moments. 

In  the  table  there  are  shown  the  results  of  the  computations  of  the  value^ 
for  15  centrifugal  atomizers,  the  geometrical  characteristics  of  which  changed  fi'ora 
2.09  to  7.50  and  the  drop  in  the  pressure  on  them  6  •  lo5— 30  •  n/ra^, 
FTD-TT-65-176/l^2  7 


From  the  table  it  la  seen  that  the  value  jj  practically  does  not  depend  on  the 
geometrical  cbaracterlatlcs  of  the  atomizers  nor  on  the  system  of  their  operation, 
and  it  varies  with  the  maximum  scattering  by  not  more  than  lOJt  around  the  average 
value  =  0.190. 


Table 


Original  Data  and  Results  of 


Determining  Parameter 


.\t  onuTa 

reoMCTpimecKaa 

xapaKTcpiicTuxa 

^lopcyNKH 

nepensA  Aaiuicuiiii 
Ha  4>^pcyiiKe  Apj,- 
•  10“*.H/ai> 

Wiicno  OTo6paHHux 
Kane;ib 

P 

1. 

4,72 

'  9 

2835 

0, 198 

2 

3,00 

10 

4521 

0|l71 

3 

G,31 

13,5 

3183- 

0, 182 

4 

5,45 

7.6 

1543 

0.186 

5 

3,65 

20 

3713 

0,206  ‘ 

6 

4,61 

30 

2603 

0,188 

7 

3,65 

10 

4008 

0.192 

8 

5, 13 

12 

4996 

oi  193 

9 

7  so 

14 

8438 

0,199 

10 

3,66 

28 

5404 

0, 192 

]  I 

3,92 

10 

3448 

0, 192 

12 

4.25 

6 

6244 

0^  186 

13 

3,45 

to 

5747 

0.168 

U 

16 

2,09 

4,34 

12 

9 

8246 

3384 

0,'  180 

0,205 

1)  ordinal  number  of  experiments;  2)  geometrical  charac¬ 
teristics  of  the  atomizer  A;  3)  drop  in  pressuire  on  the 
atomizer^  p^  •  10-‘,  n/m^;  4}  numoer  of  selected  drops 
"o.  ^ 


After  deteimlning  the  numerical  value  of  the  parameter  it  was  clear  how  well 
the  selected  function  of  distribution  describes  the  existing  data  from  observations. 
In  Figure  1  there  is  shown  for  a  single  atomizer  the  emperical  distribution  D  in  the 
form  of  a  graduated  curve  (histogram)  and  the  theoretical  distribution  (6)  with  the 
computation  for  this  atomizer  of  the  value In  Figure  2  there  is  given  for  this 
same  atomizer  the  composition  of  the  integral  curves  of  the  dlstrlbutioni  experi¬ 
mental  F(D)  and  theoretical  ^(D).  As  is  known,  F(D)  =  n(D)/no  and 


0(D)  =  <p(D)<iD=  ^D,n„exp(-pD^„/D) 


0 

(1  +?)exp(— 1 

3) 

pD„,,/D)  y 

[-+  '  ■ 

n  1  ■  —  UL/  — 

jD»(l  +  p)exp(-p) 
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Figure  1.  Histogram  ^^(0^)  and  graph 
of  function  corresponding  to  it  of  dis¬ 
tribution  density  ^(D)  for  one  centri¬ 
fugal  atomizer. 


In  Figure  3  there  is  presented  the 
griph  the  generalized  function  of  the  den¬ 
sity  of  the  distribution  of  the  dr^pj  by 
ohe  dimensions,  i.  e.,  the  graph  con¬ 
structed  in  accordance  with  the  equation 
(2)  with  the  values  =  0.190, 

(subscript  cp  =  average)  and  the  points 


i 

rT" 

» 

r~ 

I  j 

I 

1 

1 

io  i>o  60  so  no  no  o 


Figure  2.  Graph  of  tlje  experimental 
F( D)  and  theoretical  ^  ( D)  of  the  in¬ 
tegral  Auctions  of  the  distribution 
of  the  drops  of  one  centrifugal  atom¬ 
izer. 
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Figure  3.  Graph  of  the  general  ..  ; 
function  of  the  rilut  Lon  of 
hroos  in  accordan'^c  with  sizes 
(  j  )  *or  centrifugal  atomizers. 


are  plotted  according  to  the  experimental  data  for  atomizers  which  stand  out  most  by 
their  gecmietrical  characteristics  in  the  systems  of  operation  (Experiments  6,9, '4, 12) 
The  data  of  the  table  and  the  graphs  (Figures  1-3)  are  evidence  of  the  pos¬ 
sibility  of  using  equation  (?)  for  computing  the  spectra  of  atomization  by  centrifu- 
gal  atomizers  with  the  value  of  the  parameter ^  equal  to  0.190. 

For  determining  the  numerical  value  ir.  the  case  of  other  qualitatively  dif¬ 

ferent  methods  of  atomization,  it  is  necessary  to  perform  the  respective  experimei  t. 
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The  relationships  (5)  and  (?)  make  it  possible  to  compute  the  spectra  of 
atomization  by  centrifugal  atomizers  in  accordance  with  their  geometrical  dimensions, 
the  given  system  of  operation,  and  the  physical  properties  of  the  liquids,  and  can 
be  used  in  different  fields  of  technology. 

Symbols 

^  =  density 

=  dynamic  viscosity 
6  =  surface  tension 

V  =  speed  of  flow  of  liquid  from  the  atomizer 
D  =  diameter  of  the  drop 
P  diameter  of  drop 

Uq  =  overall  number 

dn  =  number  of  drops  of  the  size  from  o  to  0  +  do  or  from  D  to  D  +  dD 
^  experimentally  determined  parameter  of  the  function  of  dlstrioution 
S  =  entropy 

po=  drop  in  pressure  on  the  atomizer 

-  thlcloiess  of  the  foam  on  the  liquid  at  the  section  of  the  atomizer  nozzle 
nj  =  number  of  drops  fallen  into  the  ^th  interval  ith  interval  of  grouping 

m  =  number  of  Intervals  of  grouping 

n(D)  =  number  of  drops  the  dimension  of  which  is  less  than  D 
The  subscripts  *,  and  f  represent  respectively,  liquid  and  gas. 
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Summary 


Proceeding  from  the  Second  Law  of  Thermodynamics,  the  validity  of  Tresh'c 
assumption  (5)  for  the  derivation  of  the  function  (2)  of  droplet  size  distribution 
is  shown.  The  mathematical  statistic  methods  are  used  for  obtaining  the  analytical 
relations  between  the  parameters  of  this  function.  This  allowed  very  simple  predict¬ 
ing  formula  (9)  and  (I3)  for  the  deteraination  of  the  parameter  spectra  of  the  pul¬ 
verized  liquid.  From  the  treatment  of  the  experimental  data  the  value  of^=  0.19 
for  centrifugal  pulverizers  was  found. 

Figures  2  and  3  show  how  distribution  functions  (2)  and  (I4)  agree  with  the  ob¬ 
servation  data. 
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